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ABSTRACT 

We use a combination of high resolution TV-body simulations and semi-analytic tech- 
niques to follow the formation, the evolution and the chemical enrichment of galaxies 
in a ACDM Universe. We model the transport of metals between the stars, the cold 
gas in galaxies, the hot gas in dark matter haloes, and the intergalactic gas outside 
virialized haloes. We have compared three different feedback schemes. The 'retention' 
model assumes that material reheated by supernova explosions is able to leave the 
galaxy, but not the dark matter halo. The 'ejection' model assumes that this material 
leaves the halo and is then re-incorporated when structure collapses on larger scales. 
The 'wind' model uses prescriptions that are motivated by observations of local star- 
burst galaxies. We require that our models reproduce the cluster galaxy luminosity 
function measured from the 2dF survey, the relations between stellar mass, gas mass 
and metallicity inferred from new SDSS data, and the observed amount of metals in 
the ICM. With suitable adjustment of the free parameters in the model, a reasonable 
fit to the observational results at redshift zero can be obtained for all three feedback 
schemes. All three predict that the chemical enrichment of the ICM occurs at high red- 
shift: 60-80 per cent of the metals currently in the ICM were ejected at redshifts larger 
than 1, 35-60 per cent at redshifts larger than 2 and 20-45 per cent at redshifts larger 
than 3. Massive galaxies are important contributors to the chemical pollution: about 
half of the metals today present in the ICM were ejected by galaxies with baryonic 
masses larger than 10 10 h~ Y M@. The observed decline in baryon fraction from rich 
clusters to galaxy groups is reproduced only in an 'extreme' ejection scheme, where 
material ejected from dark matter haloes is re-incorporated on a timescale compara- 
ble to the age of the Universe. Finally, we explore how the metal abundance in the 
intergalactic medium as a function of redshift can constraint how and when galaxies 
ejected their metals. 

Key words: galaxies: formation - galaxies: evolution - galaxies: intergalactic medium 
- galaxies: stellar content - galaxies: cluster: general 



1 INTRODUCTION 

A^-body simulations have shown that the baryon fraction 
in a rich cluster does not change appreciably during its 
evolution JWhite et al.lll993lh Clusters of galaxies can thus 
be considered as closed systems, retaining all information 
abou t their past st ar formation and metal production histo- 
ries jRenzinilll99'it) . This suggests that direct observations 
of elemental abundances in the intra-cluster medium (ICM) 
can constrain the history of star formation in clusters, the 
efficiency with which gas was converted into stars, the rel- 
ative importance of different types of supernovae, and the 



* Email: gdelucia@mpa-garching.mpg.de 



mechanisms responsible for the ejection and the transport 
of metals. 

The last decade has witnessed the accumulation 
of a large amount of data on the chemical compo- 
sition of the intra-cluster gas jMush otzkv et alJ Il996l : 
iDe Grandi fc MolendfcOOltlEttori et al.l2002j) . X-rav satel- 
lites have provided a wealth of information about the abun- 
dances of many different elements. These studies have shown 
that the intra-cluster gas cannot be entirely of primordial 
origin - a significant fraction of this gas must have been pro- 
cessed in the cluster galaxies and then transported from the 
galaxies into the ICM. 

The total amount of iron dispersed in the ICM is of 
the same ord er of magnitude as the mass of iron locked in 
the galaxies iRenzini et al.lll993f) . Observational data sug- 
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gest that the mean metallicity of the ICM is about 0.2-0.3 
solar bot h for nearby (|Edge_Jg_Ste^arj|l99l|) and for distant 
clusters ijMushotzkv fc LoeweTst^!nl^97l v 

Various physical mechanisms can provide viable ex- 
planations for the transfer of metals from the galax- 
ies into the ICM, for example ejection of enric hed ma- 
terial from mergers of proto-galactic fragments (iG iifidinl 
119981); tidal/ram pre ssure stripping (iGurm JgGott I Il972t 



.Fukumoto fc Ikeuchilll99f)l: iMori fc Burkertll2000Fl: galactic 
outflowsllLaraor i£^Dmer£^^ 

Il997l: IWiebe. Shustov fc Tutukovlll999l) . The actual contri- 
bution from each of these mechanisms is still a matter of 
debate. Observational data suggest that metals are most ef- 
fectively transferred into t he ICM by inte rnal mechanisms 
rather than external forces. lRenzini| (|l997l ) has argued that 
ram pressure cannot play a dominant role, because this pro- 
cess would operate more efficiently in high velocity disper- 
sion clusters. A correlation between the richness of the clus- 
ter and its metal content is not supported by observations. 

In recent years, supernova-driven outflows have re- 
ceived increasing attention as the most plausible explana- 
tion for the presence of me tals i n the ICM. It was originally 
suggested bv lLarsonl lll974t) and lLarson fc Dinersteinl jl975t) 
that the fraction of mass and hence of metals driven from 
a galaxy increases with decreasing galactic mass, because 
lower mass galaxies have shallower potential wells. This nat- 
urally establishes a metallicity-mass relationship that is in 
qualitative agreement with the observations. It also predicts 
the chemical pollution of the ICM as a side-effect of the 
outflow. 

This outflow scenario is not without its problems, 
however. It has been shown in a number of papers 
llDavid. Forman fc Jonesl Il99ll : iMatteucci fc Gibson! 1199,4 



Gibson fc Matteuccil fl997t iMoretti. Portinari fc Chiosil 
2003) that if a standard IMF and chemical yield is assumed, 
it is very difficult to account for the total amount of metals 
observed in rich clusters. Some authors have suggested 
that cluster ellipticals may form with a non-standard 
'top-heavy' IMF. This would alleviate the metal budget 
problems and also explain the 'tilt' of the fundamental 
plane, i.e. the increas e in galaxy mass-to-light ratio with 
increasing luminosity (IZepf fc Silklll996l ; IChiosi et al1ll998t 
iPadmanabhan et all 12003) ! From a theoretical point of 
view, a skewness of the IMF towards more massive stars at 
higher redshift might be expected from simple arguments 
related to the Jeans scale and to the scale of mag netic 
support against gravitational collapse jLarsonl Il998l) . On 
the other hand, many authors have argued that there is 
very little real observational evidence that the IMF does 
vary, at least between different regions of our own Galaxy 
jHernandez fc FerraralEool: iMassevI 11991: iKrouna fc Boill 
2002). In this analysis, we will simply sweep the IMF issue 
under the carpet by treating the chemical yield as one of 
the parameters in our model. 

Many interesting clues about metal enrichment at high 
redshift have been found by studying Lyman break galax- 
ies (LBGs). Detailed studies at both optical and infrared 
wavelengths have shown that at redshifts ~ 3 , the metal- 
licity of LBGs is re l ativel y high (in the range 0.1-0.5 Zq) 
l|Pettini et alJ 1200(1 1200 ll) . Studies of the spectral energy 
distribution of these objects have shown, somewhat surpris- 
ingly, that 20 per cent of these galaxies have been forming 



stars for more than 1 Gyr iShaplev et al .1200 ll) . This pushes 
the onset of star formation in these objects to redshifts in 
excess of 5. Although there seems to be a general consen- 
sus that star formation activity (and hence the chemical 
pollution of the interstellar medium) must have started at 
high redshift, the question of which galaxies are responsi- 
ble for this pollution is still controversial. Some theoretical 
studies show that elliptical galaxies must have played an 
important role in establishing the observed abundance of 
the ICM, but these studies often require, as we noted be- 
fore, an initial mass function (IMF) that is skewed towards 
more massive stars at high redshift. Most such modelling 
has also not b een carried out in a fully cosm ological context 
(see, however. lKauffmann fc Charlod l|1998l) f or an approach 
closer to that of this paper). Other studies jGarnetdl2002l) 
suggest that dwarf galaxies have been the main contribu- 
tors to the chemical pollution of the inter-galactic medium 
(IGM). 

In this paper we use a combination of high-resolution 
./V-body simulations of the formation of clusters in a ACDM 
Universe and semi-analytic techniques to follow the enrich- 
ment history both of galaxies and of the ICM. We test that 
our model is able to reproduce observations of the number 
density, the stellar populations and the chemical properties 
of cluster galaxies, as well as the metal content of the ICM. 
We then study which galaxies were primarily responsible for 
polluting the ICM and when this occurred. 

The paper is structured as follows: in SecElwe describe 
the simulations used in this work; in Sec. [3] we summarise 
the semi-analytic technique we employ; while in Sec. 0|we 
give a detailed description of the prescriptions adopted to 
parametrise the physical processes included in our model. 
In Sec. we describe how we set the free parameters of our 
model and we show the main observational properties that 
can be fit. Sec.HJand Sec. |7|present the main results of our 
investigation on the chemical enrichment history of the ICM 
and the IGM, and investigate two observational tests that 
may help to distinguish between different feedback schemes. 
Our conclusions are presented in Sec. |H1 



2 iV-BODY SIMULATIONS 

In this study we use a collisionless simulation of a clus- 
ter of galaxies, generated using the 'zoom' technique 
llTormen. Bouchetfc Whit.fi 119971 : Isee also Katz fc White! 
1993). As a first step, a suitable target cluster is selected 
from a previously generated cosmological simulation. The 
particles in the target cluster and its immediate surround- 
ings are traced back to their Lagrangian region and replaced 
with a larger number of lower mass particles. These par- 
ticles are then perturbed using the same fluctuation field 
as in the parent simulation, but now extended to smaller 
scales (reflecting the increase in resolution). Outside the 
high-resolution region, particles of variable mass, increasing 
with distance, are displaced on a spherical grid whose spac- 
ing grows with distance from the high-resolution region and 
that extends to the box size of the parent simulation. This 
method allows us to concentrate the computational effort on 
the cluster of interest and, at the same time, to maintain a 
faithful representation of the large-scale density and veloc- 
ity of the parent simulation. In the following, we will refer 
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to haloes that include low-resolution particles as 'contami- 
nated' haloes. We will exclude these haloes in the analysis. 

The cluster simulation used in this work was carried 
out by Ba rbara Lanzoni a s part of h er PhD thesis and is de - 
scribed in lLanzoni et alJ i2003l) and lDe Lucia et alJ |2003). 
We also use a re-simulation of a 'typical' region of the Uni- 
verse, carried out by Felix Stoehr as part of his PhD thesis. 
This used the same re-simulation technique as the cluster 
simulation. In both, the fraction of haloes contaminated by 
the presence of low-resolution particles is ~ 3 per cent, all 
near the boundary of the high resolution region. 

The parent simulation employed is, in both 
cases, the Very Large Simulat ion (VLS) c arried 
out by the Virgo Consortium (fienkins et alJ 1200 it 
Isee also Yoshida. Sheth fc Diaferiol 1200 J). The simulation 
was p erformed using a parallel P3M code iMacfarland et aD 
1998) and followed 512 3 particles with a particle mass of 
7 x 10 10 h- 1 M© in a comoving box of size 479/i _1 Mpc 
on a side. The parent cosmological simulation is charac- 
terised by the following parameters: Qq = 0.3, = 0.7, 
spectral shape F = 0.21, h = 0.7 (we adopt the convention 
Ho = 100 h kms -1 Mpc^ 1 ) and spectral normalisation 
as = 0.9. 

The numerical parameters of the simulations used in 
this work are summarised in Table Q 



ies were previously central galaxies of another halo, which 
then merged to form the larger object. Because the core of 
the parent halo is still intact, the positions and velocities of 
these halo galaxies can be accurately determined. Note that 
gas is no longer able to cool onto halo galaxies. 

Dark matter subhaloes lose mass and are eventually de- 
stroyed as a result of tidal stripping effects. A galaxy that is 
no longer identified with a subhalo is called a satellite. The 
position of the satellite is tracked using the position of the 
most bound particle of the subhalo before it was disrupted. 
Note that if two or more subhaloes merge, the halo galaxy of 
the smaller subhalo will become a 'satellite' of the remnant 
sub halo. 

ISpringel et al] i200 J) show that the inclusion of sub- 
haloes results in a significant improvement in the cluster lu- 
minosity function over previous semi-analytic schemes, and 
in a morphology-radius relationship that is in remarkably 
good agreement with the observational data. This improve- 
ment is mainly attributed to a more realistic estimate of the 
merger rate: in the standard scheme too many bright galax- 
ies merge with the central galaxy on short time-scales. This 
produces first-ranked galaxies that are too bright when com- 
pared with observational data. It also depletes the cluster 
luminosity function around the 'knee' at the characteristic 
luminosity. 



3 TRACKING GALAXIES IN A-BODY 
SIMULATIONS 

The prescriptions adopted for the different physical pro- 
cesses included in our model are described in more detail 
in the next section. In this section we summarise how the 
semi-analytic model is grafted onto the high resolution N- 
body simulation. The te chniques we employ in this work are 
similar to those used bv lSpringel et al] ( 2001). 

In standard semi-analytic models, all galaxies are lo- 
cated within dark matter haloes. Haloes are usually iden- 
tified in a simulation using a standard friends-of-friends 
(FOF) algorithm with a linking length of 0.2 in units of 
the mean particle separation. The novelty of the analysis 
technique developed by Springel et al., is that substructure 
is also tracked within each halo. This means that the dark 
matter halo within which a galaxy forms, is still followed 
even after it is accreted by a larger object. The algorithm 
used t o identify subhaloes (SUBFIND) is described in de- 
tail bv lSpringel et all (l200lf) . The algorithm decomposes a 
given halo into a set of disjoint and self-bound subhaloes, 
identified as locally o verdense regions in the density field of 
the background halo. iDe Lucia et al] i2003h have presented 
an extensive analysis of the properties of the subhalo pop- 
ulation present in a lar ge sample of ha l oes w ith a range of 
different masses. As in IDe Lucia et al] B2003I1 . we consider 
all substructures detected by the SUBFIND algorithm with 
at least 10 self-bound particles, to be genuine subhaloes. 

An important change due to the inclusion of subhaloes, 
is a new nomenclature for the different kinds of galaxies 
present in the simulation. The FOF group hosts the 'central 
galaxy'; this galaxy is located at the position of the most 
bound particle in the halo. This galaxy is fed by gas cooling 
from the surrounding hot halo medium. All other galaxies 
attached to subhaloes are called 'halo galaxies'. These galax- 



4 THE PHYSICAL PROCESSES GOVERNING 
GALAXY EVOLUTION 

Our treatment of the physical processe s driving galaxy 
evolut ion is similar to the one ad opted in iKauffmann et all 
l|l999fl and ISpringel et al.l (l200lTl . Many prescriptions have 
been modified in order to properly take into account 
the exchange of metals between the different phases. 
We also include metallicity-dependent cooling rates and 
luminosities. The details of our implementation are de- 
scribed below. The reader is referred to previous papers 
for more genera l inform a tion on semi-analytic techniques 
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4.1 Gas cooling 

Gas cooling is treat ed as in IKauffmann et aD lll999l) and 
ISpringel et al.l ]200l]) . It is assumed that the hot gas within 
dark matter haloes initially follows the dark matter distri- 
bution. The cooling radius is defined as the radius for which 
the local cooling time is equal to age of the Universe at that 
epoch. 

At early times and for low-mass haloes, the cooling ra- 
dius can be larger than the virial radius. It is then assumed 
that the hot gas condenses out on a halo dynamical time. 
If the cooling radius lies within the virial radius, the gas 
is assumed to cool quasi-statically and the cooling rate is 
modelled by a simple inflow equation. 

Note that the cooling rates are strongly dependent 
on the temperature of the gas and on its metallicity. 
We model these dep endences using the collision al ionisa- 
tion cooling curves of ISutherland fc Dopital Jl993|) . At high 

10 s K) temperatures, the cooling is dominated by the 
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Table 1. Numerical parameters for the simulations used. Both the simulations were carried out assuming a ACDM cosmology with 
cosmological parameters Qo = 0.3, Ha = 0.7, T = 0.21, erg = 0.9, and h = 0.7. In the table, we give the particle mass m p in the high 
resolution region, the starting redshift z s t a rt of the simulation, and the gravitational softening e in the high-resolution region. 



Name 


Description 


m p [H^Mq] 


■^start 




gl 


10 15 fl _1 M cluster 


2.0 x 10 9 


60 


5.0 


M2 


field simulation 


9.5 x 10 s 


70 


3.0 



bremsstrahlung continuum. At lower temperatures line cool- 
ing from heavy elements dominates (mainly iron in the 10 6 - 
10 7 K regime, with oxygen significant at lower tempera- 
tures). The net effect of using metallicity-dependent cooling 
rates is an overall increase of the brightness of galaxies, be- 
cause cooling is more efficient. This effect is strongest in low 
mass haloes. 

As noted since the work of lWhite fc Frenkl(ll99ll) . these 
prescriptions produce central cluster galaxies that are too 
massive and too luminous to be consistent with observations. 
This is a manifestation of the 'cooling flow' problem, the fact 
that the central gas in clusters does not appear to be cool- 
ing desp ite the short estimated cooling time. As in p revious 
models IKauffman n et al ] |l999l: Springcl et al. 2001) we fix 
this problem ad hoc by assuming that the gas does not cool 
in haloes with V v i T > V c u t- In our model Kut = 350 km s -1 . 

Note that following ISpringel et alJ <l200lf) . we define the 
virial radius R v i r of a FOF-halo as the radius of the sphere 
centred on its most-bound particle which has an overden- 
sity 200 with respect to the critical density. We take the 
enclosed mass M vir = 100H 2 Rl iT /G as the virial mass, and 
we define the virial velocity as V v 2 ir = GM v i T /R v i T . The mass 
of a subhalo, on the other hand, is defined in terms of the 
total number of particles it contains. The virial velocity of a 
subhalo is fixed at the velocity that it had just before infall. 

Although our cooling model is extremely simplified, it 
has been shown that it produces results that are in good 
agreement with more detailed iV-body + hydrodynami- 
cal simulations that ad opt the same physics jYoshida et alJ 
l2002tlHellv et alJl2003ft . 



4.2 Star formation 

The star formation 'recipes' that are implemented in semi- 
analytic models are always subj ect t o considerable uncer - 
tainty. In lKauffmann et all <ll999l) and lSpringel et al] (l200lfl 
it is assumed that star formation occurs with a rate given 
by: 



1p = aMcold/tdyr, 



(1) 



where M co id and tdyn = -Rvir/lOKdr are the cold gas mass 
and the dynamical time of the galaxy respectively, and a 
represents the efficiency of the conversion of gas into stars. 

Previous semi-analytic models have assumed that the 
efficiency parameter a is a constant, independent of galaxy 
mass and redshift. There are observational indications, how- 
ever, that low mass galaxies conv ert gas into stars less ef - 
ficiently than high mass galaxies (Kauff mann et all 12003). 
An effect in this direction is also expected in detailed mod- 
els of the effects of supernovae feedback on the inte rstellar 
medium jMcKee fc Ostrikerl fl977l: lEfstathioul l2000l) . More 
importantly, perhaps, the above prescription leads to gas 



fractions that are essentially independent of the mass of 
the galaxy. In practice, we know that gas fractions increase 
from ~ 0.1 for luminous spirals like our own Milky Way, to 
more than 0.8 for low-mass irregular galaxies (Boissi er et al] 
1200 ill . 

In this work, we assume that a depends on the circular 
velocity of the parent galaxy as follows: 

Kir 



cvo ■ 



220 km s- 



and we treat an and n as free parameters. Note also that R V [ T 
decreases with redshift for a galaxy halo with fixed circular 
velocity. This means that a galaxy of circular velocity V c 
will be smaller at higher redshifts and as a result, the star 
formation efficiency will be higher. 

4.3 Feedback 



Previous work (IKauffmann fc Chariot] Il998l : 

ISomerville fc Primackl Il999t ICole et al] 1200(7) has shown 
that feedback processes are required to fit the faint end 
slope of the luminosity function and to fit the observed 
slope of the colour-magnitude relation of elliptical galaxies. 
The theoretical and observational understanding of how the 
feedback process operates is far from complete. 

In many models it has been assumed that the feedback 
energy released in the star formation process is able to re- 
heat some of the cold gas. The amount of reheated mass is 
computed using energy conservation arguments and is given 
by: 

AAirohcatcd = ^SN SN ^j^^ ^ 

where t/sn is the number of supernovae expected per solar 
mass of s tars fo rmed (6.3 x 10" 3 Mq 1 assuming a universal 
Salpete 3 (I1955T) IMF) and £ S n is the energy released by 
each supernova (~ 10 51 erg). The dimensionless parameter 
e quantifies the efficiency of the process and is treated as 
a free parameter. Note that changing the IMF would also 
change the amount of energy available for reheating the gas. 

One major uncertainty is whether the reheated gas 
leaves the halo. This will depend on a number of factors, 
including the velocity to which the gas is accelerated, the 
amount of intervening gas, the fraction of energy lost by ra- 
diative processes, and the depth of the potential well of the 
halo. 

On the observational side, evidence in support of the 
existence of ou tflows from galaxies has grown rapidly in 
the last year s llHeckman et aljll995t iMarlowe et al.l Tl995l : 
iMart id 119961. In many cases, the observed gas velocities 
exceed the escape velocity of the parent galaxies; this ma- 
terial will then escape from the galaxies and will be in- 
jected into the intergalactic medium (IGM). Observations of 
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galactic-scale outflow s of gas i n active star-forming galaxies 
jLehnert fe HeckmarJ |l996|: |Dahlem 1 Weaver fc Heckmanl 
ll998l:lHeckman et all200nHrIeckmarJ2002F l suggest that out- 
flows of multiphase material are ubiquitous in galaxies in 
which the global star-formation rate per unit area exceeds 
roughly 1O -1 M0 yr _1 kpc -2 . Different methods to estimate 
the outflow rate suggest that it is comparable to the star for- 
mation rate. The estimated outflow speeds vary in the range 
400-800 km s" 1 and are independent of the galaxy circular 
velocity. These observational results indicate that the out- 
flows preferentially occur in smaller galaxies. This provides a 
natural explanation for the observed relation between galaxy 
luminosity and metallicity. 

An accurate implementation of the feedback process is 
beyond the capabilities of present numerical codes. As a re- 
sult, published simulation results offer little indication of 
appropriate recipes for treating galactic winds. In this pa- 
per we experiment with three different simplified prescrip- 
tions for feedback and study whether they lead to different 
observational signatures: 

• in the retention model , we us e the prescriptions adopted 
by iKauffmann fc Charlod dl998tl and assume that the re- 
heated material, computed according to Eq. [5] is shock 
heated to the virial temperature of the dark halo and is 
put directly in the hot phase, where it is then once more 
available for cooling. 

• In the ejection model, we assume that the material re- 
heated by supernovae explosions in central galaxies always 
leaves the halo, but can be later re-incorporated. The time- 
scale to re-incorporate the gas is related to the dynamical 
time-scale of the halo by the following equation: 

AM back = 7 • Mejeoted ' ■ At (3) 

where AM Dac k is the amount of gas that is re-incorporated 
in the time-interval At; M c j G ctcd is the amount of mate- 
rial in the ejected component; ii v j r and VW are the virial 
radius and the virial velocity of the halo at the time the 
re-incorporation occurs; 7 is a free parameter that controls 
how rapidly the ejected material is re-incorporated. Note 
that the material ejected in this way is not available for 
cooling until it is re-incorporated in the hot component. 

For all the other galaxies (halo and satellite galaxies), we 
assume that the material reheated to the virial temperature 
of the subhalo is then kinematically stripped and added to 
the hot component of the main halo. 

• In the wind model, we adopt prescriptions that are mo- 
tivated by the observational results. We assume that only 
central galaxies residing in haloes with a virial velocity less 
that Vcrit can eject outside the halo. The outflow rates from 
these galaxies are assumed to be proportional to their star 
formation rates, namely: 

Mw = c • V (4) 

Observational studies give values for c in the range 1-5 
jMartinlll999h while values in the range 100-300 km s are 
reasonable for Vcrit <Heckmarj 2002,. We treat Vcrit and c as 
free parameters. We also assume that the ejected material 
is re-incorporated as in the ejection scheme. If the condi- 
tions for an outflow are not satisfied, the reheated material 
(computed according to Eq.|2) is treated in the same way as 
in the retention model. The wind model is thus intermedi- 



ate between the retention and the ejection schemes. Satellite 
galaxies are treated as in the ejection scheme. 

4.4 Galaxy mergers 

In hierarchical models of galaxy formation, galaxies and 
their associated dark matter haloes form through merging 
and accretion. 

In our high resolution simulations, mergers between 
subhaloes are followed explicitly. Once a galaxy is stripped of 
its dark halo, merging timescale s are estimated using a sim - 
ple dynamical friction formula: feinnev fc Tremainel ll987'l 

_ 1 /(e) y vir ift r 

friction - 2 c GMsatlnA 

iNavarro. Frenk fc White! Il995l> show that this analytic 
estimate is a good fit to the results of numerical simulations. 
The formula applies to satellites of mass M sa t orbiting at a 
radius _R v i r in a halo of virial velocity VVir- f(e) expresses the 
dependence of the decay on the eccent ricity of the orbit an d 
is well approximated by /(e) ~ e 0,78 llLacev fc Coldll993ll : 
C is a constant ~ 0.43 and InA is the Coulomb logarithm. 
We adopt the average value < /(e) >~ 0.5, computed by 
iTormenl |l997), and approximate the Coulomb logarithm 
with InA = (1 + M v i r /M sa t). For the satellite galaxy mass 
we use the value of M V [ T corresponding at the last time the 
galaxy was a central galaxy (either of a halo or of a subhalo). 

When a small satellite merges with the central galaxy, 
its stellar mass and cold mass are simply transferred to the 
central galaxy and the photometric properties are updated 
accordingly. In particular we transfer the stellar mass of the 
merged galaxy to the bulge of the central galaxy and up- 
date the photometric properties of this galaxy. If the mass 
ratio between the stellar component of the merging galaxies 
is larger than 0.3, we assume that the merger completely de- 
stroys the disk of the central galaxy producing a spheroidal 
component. In addition we assume that the merger con- 
sumes all the gas left in the two merging galaxies in a single 
burst. The stars formed in this burst are also added to the 
bulge. Note that since the galaxy is fed by a cooling flow, it 
can grow a new disk later on. 

4.5 Spectro— photometric evolution 

The photometric properti es of our model galaxies ar e calcu- 
lated using the models of iBruzual fc Charlod lll993l) . which 
include the effect of metallicity on the predicted luminosi- 
ties and colours of a galaxy. The stellar population synthesis 
models are used to generate look-up tables of the luminos- 
ity of a single burst of fixed mass, as a function of the age 
of the stellar population and as a function of its metallic- 
ity. When updating the photometric properties of our model 
galaxies, we interpolate between these tables using a linear 
interpolation in t and log Z. It is assumed that stars form 
with the same m etallicity as the cold gas. We have adopted 
a lSalpeterl lll955l) IMF with upper and lower mass cut-offs 

of 1 00 and 0.1 Mp. 

ICharlot. Worthev fc Bressanl (I1996T) have demonstrated 
that for a given IMF and star formation history, the broad- 
band colours produced by different stellar population codes 
differ by only a few tenths of a magnitude. The most im- 
portant sources of uncertainty in our model predictions are 
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thus the IMF (which affects the luminosity quite strongly) 
and the associated yield (which influences the colours). 



4.6 Dust extinction 

Attenuation of starlight by dust affects the colours of galax- 
ies. The properties of dust are dependent on a number of 
factors such as the star formation rate, that regulates the 
rate of creation, heating and destruction of dust grains and 
the distribution of dust and metals within gas clouds. For a 
single galaxy, all these factors can be taken into account, and 
it is then possible to model the effect of dust on the galaxy's 
spectrum. However, the level of detail that is required goes 
far beyond the capabilities of our present code. We there- 
fore adopt a dust model that is based on the macroscopic 
properties of galaxies, i.e. luminosity and inclination. 

This model has been used in previous work 
jKauffmann et alJ llQflfll ISomerville fc Primackl Il99flft and 
is bas ed on the observational results bv lWang fc Heckmanl 
( 1996), who studied the correlation between the face-on op- 
tical depth of dust in galactic discs and the total luminosity 
of the galaxy. IWang fc Heckmanl il99rJl find that this can 
be expressed as: 



TB = TB,: 



Lt 



where Lb is the intrinsic (unextincted) blue luminosity and 
Lb,* is the fiducial observed blue lumin osity of a Schechter 
L„ ga laxy (M*(B) = -19.6 + 51og/i). IWang fc Heckmanl 
(1996) find that a relation with tb,» = 0.8 ±0.3 and f3 ~ 0.5 
fit their data very well. 

We use tb,» = 0.8, /3 = 0.5 and relate the B-band op- 
ti cal depth to the other bands usin g the extinction curve 
of lCardelli. Clayton fc Mathisl <ll989F) . We also assign a ran- 
dom inclination to each galaxy and apply the dust correction 
only to its disc component (i.e. we assume that the bulge is 
not affected by dust) using a 'slab' geometry: 



A\ = -2.51oi 



1 - e~^ se 
T\sec9 



where 9 is the angle of inclination to the line of sight. All 
the magnitudes and the colours plotted in this paper include 
the effects of dust extinction, unless stated otherwise. 



4.7 Metal routes 

We assume that a yield Y of heavy elements is produced per 
solar mass of gas converted into stars. All the metals are 
instantaneously returned to the cold phase (note that this 
means that we are assuming a mixing efficiency of 100 per 
cent). We also assume that a fraction (R) of the mass in 
stars is returned to the cold gas. 

Metals are then exchanged between the different gas 
phases depending on the feedback model (see Sec. 14.31 . In 
the retention model, the metals contained in the reheated 
gas are put in the hot component and can subsequently cool 
to further enrich the cold phase. In the ejection and wind 
schemes, the metals can be ejected outside the halo. The 
ejected metals are re-incorporated into the hot halo gas on a 
halo dynamical time (see Eq.[3J. In the wind scheme, metals 
are only ejected for haloes with V v it < Vcrft • In more massive 



systems, the metals contained in the reheated gas are mixed 
with the hot component of the main halo. Some simulations 
jMac Low fc Ferraralll999l) suggest that low-mass galaxies 
may lose essentially all their metals, but it is difficult for 
these galaxies to 'blow-away' their interstellar medium. For 
simplicity, we do not consider schemes in which metals are 
selectively ejected from galaxies. 

When a satellite galaxy merges, its stars, cold gas and 
metals are simply added to those of the central galaxy. If a 
major merger occurs, all the gas is consumed in a starburst 
and the metals are ejected into the hot phase in the retention 
scheme and into the ejected phase in the other two schemes 
(in the wind scheme ejection occurs only in galaxies that 
reside in haloes with V v ii > Vcrft). 

In Fig. we sketch the routes whereby mass and 
metals are exchanged in the model. We now write down 
the equations that describe the evolution of the mass of 
the four reservoirs shown in Fig. Q All central galaxies 
have 4 different components: stars, cold gas, hot gas and 
an 'ejected' component. The equations describing these 4 
components (in the absence of accretion of external matter) 



M atars = (1 - R) ■ V 

A/hot = — Mcool + Mback + X^aat Af reheated 
Mcold = +Mcool - (1 - R) ■ lp ~ Mout 
A/ejected = + M ou t — Mback 

where M out is the rate at which the cold gas is ejected 
outside the halo and is given by Eq. for the ejection 
scheme and Eq. [I] for the wind scheme; V is the star 
formation rate given by Eq. M coo i is the cooling rate; 
A/back is the rate at which gas is reincorporated from the 
ejected material (Eq. El; J2 B s.t -^reheated is the sum of all 
material that is reheated by the satellite galaxies in the 
halo. 

Satellite galaxies have no hot gas or ejected components 
and M coo i = 0. The equations for the satellite galaxies are: 

M stars = (1 - R) ■ V 

M co ld = + M CO ol — (1 — R) ■ 4> — M rc hcatcd 

where M re heated is computed using Eq. [5] 

From the above equations, one can easily obtain the 
corresponding equations for the evolution of the metal 
content of each reservoir of material. For central galaxies, 
the equations are: 



M s f ars = +(l-R)-^-^ c , 



M,f 



- Mcool • Zhot + Mback 1 -^ejected +X] S at [-^reh 



Afcold = +M C ool • Zhot - (1 - R) ■ >0 ■ Zedd +Y-lj} — Mout • Zcold 
Mefected = "t"Mout ' ^cold — Mback ' ^ejected 

where Z cold = M c ^i d /M co id, ^ejected = M^ ected /M ejected and 
■Zhot = Mj^t/Mhot represent the metallicities in the cold 
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re-incorporation 




Figure 1. A schematic representation of how mass is exchanged 
between the different phases considered in our models in the ab- 
sence of accretion from outside. Each arrow is accompanied by a 
name indicating the physical process driving the mass exchange. 
Metals follow the same routes as the mass. 



gas phase, in the ejected component and in the hot gas 
phase respectively. 

Note that the above equations assume that material 
is being ejected outside haloes. In our 'retention' scheme, 
no material leaves the halo and the applicable equations 
are obtained by neglecting the equation for M e j ec ted and 
substituting M back and M ou t with M re heated (given by Eq.0. 

The practical implementation of th e prescriptions is de- 
scribed in detail in lSpringel et alJi200ll) . We have stored 100 
outputs of the simulation spaced in equal logarithmic inter- 
vals in redshifts from z = 20 to z = 0. For each new snap- 
shot, we estimate the merger timescales of satellite galaxies 
that have entered a given halo. The 'merger clock' for the 
other satellites is updated and a galaxy merges with the 
central object when this time has elapsed. The merger clock 
is reset if the halo containing the satellites merges with a 
larger system. 

The total amount of hot gas available for cooling in each 
halo is given by: 

M hot = / b M vir - £ [M s « + Mil + <L cd ] (5) 

i 

where the sum extends over all the galaxies in the halo and 
/b is the baryon fraction of the Universe. In this work we 
use the value / b = 0.15 as s uggested by the CMB exper- 
iment WMAP dSpereel et al.ll2003ft . When going from one 
snapshot to the next, three things may happen: 

• a new halo may form. In this case its hot gas mass is 
initialised to the value /bM,i r ; 

• the virial mass of the halo may increase because of ac- 
cretion of 'diffuse' material. In this case the accreted baryons 
are effectively added to the hot component using Eq. |S] 

• two haloes can merge. In this case both the ejected 
component and the hot component of the lower mass halo 
are added to the hot component of the remnant halo. The 
ejected component of the remnant halo remains outside the 
halo but is reduced by a factor given by Eq. [3] 

The cooling rate is assumed to be constant between 
two successive simulation outputs. The differential equations 



given above are solved using smaller time-steps (50 between 
each pair of simulation snapshots) . 

4.8 Model normalisation and influence of the 
parameters 

As in previous w o rk llKauffmann fc Charlotl Il998l : 
ISomerville fc Primackl Il999f) . the free parameters of 
the model are tuned in order to reproduce the observed 
properties of our Galaxy (see Table |3J. We also check 
that we get the right normalisation and slope for the 
Tully-Fisher relation. 

Note that the main parameters controlling the physical 
evolution of our model galaxies are: 

• the parameters determining the star formation effi- 
ciency: ao and n; 

• the parameters determining the feedback efficiency: e in 
the ejection model, and c and Krit in the wind model; 

• the amount of metals produced per solar mass of cold 
gas converted into stars: Y; 

• the gas fraction returned by evolved stars: R. 

• the parameter 7 that determines how long it takes for 
the ejected gas to be re-incorporated. 

The value of the parameter R can be directly estimated 
from stellar evolution theory for a given choice of IMF. Pop- 
ulation synthesis models show that this recycled fraction is 
roughly independent of metallicity and lies in the range 0.2- 
0.45 JCole et all200rA . Note that this value can be larger for 
top-heavy IMFs. 

Note that there are only a few free parameters in our 
model. The influence of the different parameters on the ob- 
served properties of galaxies can be summarised as follows: 

• the yield Y controls the total amount of metals in the 
stars and gas. Note that the cooling rates in lower mass 
haloes are strongly dependent on the metallicity. Increasing 
the yield results in an overall increase of the luminosities of 
our galaxies and a tilt in the Tully-Fisher relation towards 
a shallower slope. 

• The star formation efficiency ao has only a weak influ- 
ence on the zero-point of the Tully-Fisher relation, but it 
has an important influence on the gas fraction of galaxies. 

• The parameter n, that parametrises the dependence 
of star formation efficiency on the circular velocity of the 
galaxy, has a strong effect on the dependence of the gas 
fraction on the mass or circular velocity of the galaxy. It has 
negligible effect on all the other observational properties. 

• The feedback efficiency has a strong influence both on 
the zero-point and the slope of the Tully-Fisher relation. 
An increase in feedback efficiency results in a decrease in 
the luminosities of galaxies and a tilt in the Tully-Fisher 
towards a steeper slope. 

• The gas fraction returned by evolved stars R has only 
a marginal influence on the gas metallicity and on the lumi- 
nosities of the galaxies in our model. 

• The parameter 7 controls how long it takes for the 
ejected gas to be re-incorporated back into a dark matter 
halo. We have experimented with 7 values in the range 0.1-1. 
If 7 = 0.1, the re-incorporation time is of order the Hubble 
time. If 7 is large, the ejection model simply reverts back to 
the retention scheme. Note that if we decrease 7, gas and 
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Table 2. Free parameters values adopted for the three models. 





«0 


n 


c 


Y 


R c 




retention 


0.09 


2.2 


0.45 


0.045 


0.35 




ejection 


0.08 


2.2 


0.15 


0.040 


0.35 




wind 


0.10 


2.5 


0.35 


0.040 


0.35 5 


150 



metals remain outside dark matter haloes longer and cooling 
rates are reduced. As a result, the feedback efficiency in our 
preferred model is smaller. In the next section, we will show 
results for 7 = 0.1, which should be considered an 'extreme' 
ejection scheme, where gas and metals remain outside the 
halo for a time comparable to the age of the Universe. In 
Sec. [7] we will explore in more detail what happens if the 
re-incorporation timescale is reduced. 

A through exploration of parameter space results in the 
parameters listed in Tabled The requirement that our mod- 
els agree with a wide range of observational data (see Sec.|5J 
allows only slight changes around the values listed in Table 

m 

Note that in all three models, we have to as- 
sume a value for the yield that is larger than the con- 
ventional value. This is consistent with other analyses, 
which have shown that the observed metallicity of the 
ICM cannot be explained usin g a standard IMF with 
a standard value of the yield | ffit«on fc Matteucclll997t 
iGibson. Loewenstein fc Mushot zky[^^9%. A recen t review 
of the problem can be found in lMoretti et alJ (I20031) . In this 
work, we will leave aside this problem. Note also that we 
will not attempt to distinguish between the heavy elements 
produced by SNII and SNIa. It is known that the latter 
are the most important contributors for Fe, while the for- 
mer mainly contribute a elements. SNII and SNIa events are 
characterised by different time-scales, so that there is a lag 
between the ejection of these elements into the interstellar 
medium. This will be studied in more detail in a future pa- 
per (Cora et al., in preparation). We will also not attempt 
to explore what would happen if the various free parameters 
listed in Table H were to depend on redshift. A redshift de- 
pendence of the star formation efficiency, for example, may 
be required to explai n the evolution in the number density 
of luminous quasars (iKauffmann fc Haehneltll200d) and to 
reproduce the o bserved properties of Lyman break galaxies 
at redshift ~ 3 iSomerville. Primack fc Faberll200il) . 

Note that the only parameter that changes significantly 
between the three different models is the feedback efficiency 
and that for the retention and wind schemes, we are obliged 
to adopt an uncomfortably high value. The value adopted 
for the parameter c also lies on the upper end of the allowed 
range. Efficient feedback is required in order to counteract 
the high cooling rates and prevent overly luminous galaxies 
from forming at the present day. If the ejected material is 
kept outside the haloes for substantial periods, then the need 
for high feedback efficiencies is not as great. Note that a 
high value of the feedback efficiency is also required in the 
wind model. This is because most of the galaxies in the 
simulation box rapidly fall below the conditions required for 
a 'wind', given the relatively low values adopted for Krit- 
Efficient feedback is then again required in order to avoid 
the formation of overly luminous galaxies. 



5 COMPARISON WITH LOCAL 
OBSERVATIONS 

In this section, we present model results for some of the basic 
observed properties of galaxies, both in and out of clusters. 
As explained in Sec. 14.81 the free parameters of our model 
are mainly tuned in order to reproduce the observed proper- 
ties of the Milky-Way and the correct normalisation of the 
Tully-Fisher relation . An extensive exploration of param- 
eter space seems to indicate that the range of acceptable 
parameters for each model is small. 

We select as 'Milky-Way type' galaxies all the objects 
in the simulation with circular velocities in the range 200- 
240 km s -1 and wit h bulge-to-disk ratios consisten t with 
Sb/Sc type galaxies iSimien fc de Vaucouleurslll986t) . More 
specifically we select all galaxies with 1.5 < AM < 2.6 
(AM = Mbuig C — Mtotai). Note that in our model we assume 
that the circular velocity of a galaxy is ~ 25 per cent larger 
than the circula r velocity of its halo. Th is is motivated by 
detailed models iMo. Mao fc Whitdll99l) for the structure 
of disk galaxies embedded in cold dark matter haloes wit h 
the universal NFW profile jNavarro. Frenk fc Whitdll997f) . 
These models show that the rotation velocity measured at 
twice the scale length of the disk is 20-30 per cent larger 
than the virial velocity of the halo. 

In our three models we find 11, 13 and 11 Milky- Way 
type galaxies (for the retention, ejection and wind model re- 
spectively) that reside in uncontaminated haloes. For these 
galaxies we obtain the stellar masses, gas masses, star for- 
mation rates and metallicities listed in Tabled These values 
are very close to what is observed for our own Galaxy, which 
has a stellar mass ~ 10 11 Mq and a total mass of cold gas 
in the range 6.5 ■ 10 9 -5.1 ■ IP 1 " M^ iKauffmann et alJll99st 
ISomerville fc Primackl Il999l) . iRocha-Pinto et alJ d2000allrl 
give an estimate for the mean SFR in the Milky Way's 
disc over the last few Gyr of the order of l-3M0yr _1 . 
This is somewhat lower than the values we find for our 
model galaxies. The Galaxy has a B-band and I-band ab- 
solute magnitude of ~ —20.5 and ~ —22.1 respectively, 
and a V-light weighted mean metallicity of ~ 0.7 sola r 
iKauffmann &: Charlod Il99£t ISomerville fc Primackl Il999^ . 
All these observed values are in reasonably good agreement 
with the values listed in Tabled] 

In Fig. [5] we show cluster luminosity functions in the 
bj-band for each of our three schemes. We find cluster lu- 
minosity functions which are in a reasonably good agree- 
ment with the comp osite luminosity f unctio n for cluster 
galaxies obtaine d by iDe Propris et all <l2003t) . A fit to a 
ISchechteil ( 1976) function gives a characteristic magnitude of 
M£ = —21.47 and a faint-end power law slope of a. = —1.22 
for the retention model. Note that resolution effects artifi- 
cially flatten the luminosity function at magnitudes > —15 
(the fit to a Schechter function is performed on the magni- 
tude interval from —16 to —22). For the ejection and wind 
schemes the results are similar. The cluster galaxies are 
slightly more luminous and the faint-end slope is somewhat 
steeper. 

Fig.[3]sliows that our model galaxies correctly reproduce 
the slope of the observed colour-magnitude relation for clus- 
ter ellip ticals. The solid line is the fit to data for Coma clus- 
ter from lBower. Lucev fc Ellis! (I1992T) . The points are model 
galaxies with an early type morphology (AM < 0.2) for our 
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Table 3. Properties of Milky-Way type galaxy in our simulation. The masses are in units of h 1 Mq and the SFR is in units of Mq yr 1 . 

M gas Mstar SFR Zstars/Z© M B -5log(h) M v - 5log(h) Af 7 -5Zog(/i) 

retention 8.03 ■ 10 9 4.67 ■ 10 10 4.57 0.80 -20.49 -21.07 -22.08 

ejection 7.53 ■ 10 9 4.94 ■ 10 10 3.76 0.82 -20.42 -21.03 -22.07 

wind 8.83 • 10 9 6.00 ■ 10 10 5.25 0.87 -20.70 -21.27 -22.30 
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Figure 2. Luminosity function in the bj-band for model galaxies 
in the retention, ejection and wind sch eme. The points re present 
the composite luminosity function by |Pe Propris et alj ( 2003 ) . 
The solid line is for the retention model, the dashed line for the 
ejection model and the dotted line for the wind model. 
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Figure 3. Colour-magnitude relation for early type galaxies in 
the simulations compared with the observationa l relat ion ob- 
tained for Coma cluster ellipticals bv lBower et al ] <1992ft . Filled 
circles are for the retention model, empty circles for the ejection 
model and crosses for the wind model. 



three different schemes. Our results confirm the conclusion 
that th e colour-magnitude i s mainly driven by metallicit y 
effects jKodama et alJ Il998l : iKauffmann fc CharlotJ Il99ct . 
Note that the scatter of the model galaxies is larger than 
the observational value: for all three mo del we find a s catter 
of ~ 0.07, while the value measured bv lBower et alJ (Il992f) 
is 0.048. 

In Fig. [I] we show the Tully-Fisher relation obtained 
for our model galaxies. We have plotted central galaxies 
that are in haloes outside the main cluster and that are 
not contaminated by low resolution particles. We apply a 
morphological cut that selects Sb/Sc-type galaxies (1.5 < 
AM < 2.6) and we select all galaxies brighter than —18 in 
the I-band. These morphological and magnitude cuts corre- 
sp ond approximate l y to t he ones defining the sample used 
bv lGiovanelli et all Il997l) . The mean observational relation 
is shown by a solid line in the figure, and the scatter about 
the relation is indicated by the dashed lines. The relation of 
Giovanelli et al. is already corrected for internal extinction. 
We therefore do not correct our I-band magnitudes for dust 
extinction in this plot. 

Note that the Tully-Fisher relations in our simulations 
have extremely small scatter. There are probably many 
sources of additional scatter that we do not account for - for 



example photometric errors and scatter in the relation be- 
tween halo circular velocity and observed line-widths. The 
slope of our predicted TF relation is also somewhat steeper 
than the observations for the retention and ejection mod- 
els (for clarity only model galaxies from the ejection scheme 
are plotted as filled circles in the figure; galaxies in the re- 
tention scheme exhibit a relation that is very similar). We 
obtain the best fit to the observed Tully-Fisher relation for 
our wind model (shown as empty circles). This is because in 
this scheme the ejected mass is systematically smaller than 
in the ejection scheme for low mass galaxies. The slope of the 
Tully-Fisher relation is strongly dependent on the adopted 
feedback prescriptions and could in principle be used to test 
different feedback models if we had better control of the 
systematic effects in converting from theoretical to observed 
quantities. 

In Fig. we compare the relation between metallicity 
and stellar mass for our simulated galaxies with the stellar 
mass-met allicity relation derived from new data from the 
Sloan Digital Sky Survey (Tremonti et al., in preparation). 
Measuring the metallicities of stars in a galaxy is difficult, 
because most stellar absorption lines in the spectrum of a 
galaxy are sensitive to both the age of the stellar popu- 
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Figure 4. Tully— Fisher relatio n for our model galaxies compared 
with the observational result bv lGiovanelli et alj Jl997f) ■ The scat- 
ter in the observational sample is shown with dashed lines while 
the thick solid line shows the 'best fit' relation given by Giovanelli 
et al. The filled circles represent model galaxies in the ejection 
model while the empty circles are used the model galaxies from 
the wind model. 



lation and its metallicity. The metallicity of the gas in a 
galaxy can be measured using strong emission lines such as 
[Oil], [OIII], [Nil], [SII] and H/3, but up to now, the avail- 
able samples have been small. Tremonti et al. have measured 
gas-phase metallicities for ~ 50, 000 emission-line galaxies 
in the SDSS and the median relation from their analysis is 
plotted as a solid line in Fig. Note that the stellar masses 
for the SDSS galaxies are measured by Tremon ti et al. using 
the same method as in lKauffmann et al] 1120031) . The dashed 
lines represent the scatter in the observed relation. The dots 
in the figure represent simulated galaxies in the retention 
model. We have only selected galaxies that reside in uncon- 
taminated haloes outside the main cluster and that have 
gas fractions of at least 10 per cent. All of these galaxies are 
star-forming and should thus be representative of an emis- 
sion line-selected sample. We obtain very similar relations 
for our two other schemes. This can be seen from the solid 
symbols in the plot, which indicate the median in bins that 
each contain ~ 400 model galaxies. Note that for the reten- 
tion and the ejection scheme, the model metallicities are sis- 
tematically lower than the observed values, although within 
la from the median observed relation. However recent work 
suggests that strong line methods (as used in Tremonti 
et al.) may systematically o verestimate oxygen abundances 
by as much as 0.2-0.5 dex llKennicutt. Bresolin fc Garnettl 
2003). Given the uncertainties in both the stellar mass and 
the metallicity measurements, the agreement between our 
models and the observational results is remarkably good. 

In Fig. we compare the metal licity ver s us ro tational 
velocity measurements published bv iGarnet u (1200211 to the 
results obtained for our models. For simplicity, we only show 



Figure 5. Metallicity-mass relation for our model galaxies 
(shown as points). The thick solid line represents the median re- 
lation between stellar mass and metallicity obtained for a sample 
of 50,000 galaxies in the SDSS (Tremonti et al., in preparation). 
The dashed lines indicate the scatter in the observed relation. The 
solid symbols represent the median obtained from our models in 
bins containing ~ 400 galaxies each. Filled circles are for the re- 
tention model, filled triangles are for the ejection model and filled 
squares are for the wind model. The error bars mark the 20th and 
80th percentiles of the distribution. 



results from the retention model. The other two schemes give 
very similar answers. Again we have selected only galax- 
ies that reside in uncontaminated haloes outside the main 
cluster and that have gas fractions of at least 10 per cent. 
The vertical dashed line in the figure shows a velocity of 
120 tons" . According to the analysis of Garnett, this marks 
a threshold below which there is a stronger dependence of 
metallicity on the potential well depth of the galaxy. Note, 
however, that the plot shown by Garnett is in linear units in 
Kot and that the 'turnover' in the mass-metallicity relation 
is much more convincing in the data of Tremonti et al. As 
can be seen, our models fit the metallicity-rotational veloc- 
ity data as well. This is not surprising, given that we obtain 
a reasonably good fit to the observed Tully-Fisher relation. 

In Fig.|7|we show a comparison between the gas fraction 
of galaxies in o ur models and the gas fractions computed by 
iGarnetu J2002). The same sample of objects as in Fig. HJ is 
plotted, both for the models and for the observations. Again 
our model agrees well with the observational data. 



6 CHEMICAL ENRICHMENT OF THE ICM 

In our simulated cluster, the metallicity in the hot gas com- 
ponent is ~ 0.26-0.30 Zp t which is in goo d agreement with 
X-ray measurements ijRenzini et al.llT993h . 

As pointed out by Renzini et al., simple metal abun- 
dances in the ICM depend not only on the total amount 
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Figure 6. Relation between metallicity and rotational velocity 
for our model galaxies (shown as p oints). The fille d squares rep- 
resent the observational data from iGarnetJ 120021) . The vertical 
dashed line in the panel corresponds to a velocity 120 km s — 1 (the 
turn-over velocity claimed by Garnett). 




Figure 7. Gas fraction as function of the B-band luminosity for 
our model galaxies (s hown as points ). Filled squares are the gas 
fractions computed by Garnett 1 2002) for the same sample of data 
shown in Fig. [S] Results are shown for the retention model, but 
are very similar for the other two schemes. 



of metals produced in stars, but also on how much dilution 
there has been from pristine gas. A quantity that is not de- 
pendent on this effect or on the total mass in dark matter in 
the cluster is the so-called iron mass-to-light ratio (IMLR) , 
defined as the ratio between the mass of iron in the ICM 
and the total B-band luminosity of cluster galaxies. 

Our simulated cluster has a total mass of 1.14 • 



Table 4. Mass-to— light ratios in the V and in the B— band for 
our three models. The units are in Mq. 

T b T v 

retention 290 230 

ejection 250 200 

wind 240 190 



10 5 h 1 Mq and a total luminosity in the B band of 
~ 10 13 L;?) . Assuming a solar iron mass fraction from 
iGrevesse. Noels fc Sauval Jl99rJ) . we find IMLR = 0.015 - 
0.020 MqLq 1 for our three models , in agreement with the 
range given bv lRenzmT et al.1 (Il993l) . 

Table|l|lists the total mass-to-light ratios of our cluster 
in the V and B-bands. The results are given for the three dif- 
ferent feedback schemes. Note that the observational deter- 
mination of a cluster mass-to-light ratio is not an easy task. 
Both the mass and the luminosity estimates are affected by 
uncertainties. Estimates based on the virial mass estima- 
tor g ive T v ~ 175-252 JCarlberg et alJll99£t I Girardi et alJ 
2000) while B-b and mass-to-light ratios are in the rang e 
200-400 h' 1 T Q llKent fc Gunnlll982l ; iGirardi et aT1l2002l) . 
It has been argued that the virial mass estimator can give 
spurious results if substructure is present in the cluster or if 
the volume sampled does not extend out to the virial radius. 
In general, estimates based on masses derived from X-ray 
data tend to give lower values. 

The good agreement between the model predictions and 
the observational results indicates that our simulation may 
provide a reasonable description of the circulation of metals 
between the different baryonic components in the Universe. 
We now use our models to generate a set of predictions for 
when the metals in the ICM were ejected and for which 
galaxies were primarily responsible for the chemical pollu- 
tion. 

Recall that we assume that metals are recycled instanta- 
neously and that the chemical pollution of the ICM happens 
through two routes: 

• in the retention scheme the reheated mass (along with 
its metals) is ejected directly into the hot component. In this 
model, the enrichment of the ICM occurs at the same time 
as the star formation. 

• In the ejection scheme the reheated material (along with 
its metals) stays for some time outside the halo and is later 
re-incorporated into the hot component. This means that 
there is a delay in the enrichment of the ICM in this scheme. 

• The wind scheme sits somewhere in between the reten- 
tion and ejection models. Metals are ejected out of the halo 
by galaxies that satisfy the outflow conditions. Otherwise, 
the metals are ejected directly into the hot gas. 

The instantaneous recycling assumption means that the 
epoch of the production of the bulk of metals in the ICM 
coincides with the epoch of the production of the bulk of 
stars. In Fig. |H| we show the average star formation rate 
(SFR) of galaxies that end-up in the cluster region and the 
average star formation rate for galaxies that end-up in the 
field (the field is defined as consisting of all haloes outside 
the main cluster that are not contaminated by low resolution 
particles). 

The SFRs in the two regions are normalised to the total 
amount of stars formed. The figure shows that star forma- 
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Figure 8. Average star formation rate as function of redshift in 
the cluster (solid line with filled squares) and in the field (dashed 
line with filled triangles). 



tion in the cluster is peaked at high redshift (~ 5) and drops 
rapidly after redshift 3. The peak of the SFR in the field oc- 
curs at lower redshift (~ 3). The decline in star formation 
from the peak to the present day is also very much shallower 
in the field. 

Integrating the SFR history of the cluster, we find that 
~ 85 per cent of the stars in the cluster formed at redshift 
larger than 2 and ~ 70 per cent at redshift larger than 3. 
The corresponding values in the field are 60 per cent and 40 
per cent. The results are similar in all the three models. 

As explained in Sec. there is no general con- 
sensus on which galaxies were responsible for enriching 
the ICM. Several theoretical studies have proposed that 
mass loss from elliptical galaxies is an effective mech- 
anism for explaining the observed amount of metals in 
the ICM. Some of these studies require an IMF that 
is skewed towards mo r e ma s sive stars at high r edshift 
jMatteucci fc VettolarJ Il988l: iGibson fc Matteucci Il997t 
iMoretti et alj I2003I) . iKauffmann fc CharlotJ dl9981) used 
semi-analytic models to model the enrichment of the ICM 
and found that a significant fraction of metals come from 
galaxies with circular velocities less than 125kms _1 . There 
have also been a number of more recent attempts to 
address this question in using hydrodynamic simulations 
jAeuirre et alJl200ltlSpringel fc Hernauisdl2003l) 

In a recent observational study. iGarnettT <|2002T) has 
analysed the dependence of the so-called 'effective yield' on 
circular velocity for a sample of irregular and spiral galaxies 
(the data are compared to our model galaxies in Fig. |BJ, 
concluding, as iLarsonl (Il974f) had earlier, that galaxies with 
Vrot ^ 100 — 150 km s -1 lose a large fraction of their SN 
ejecta, while galaxies above this limit retain most of their 



metals. As we have noted previously, this observational sam- 
ple is limited both in number and in dynamic range. 

Our model predictions are given in Fig.|U] The top panel 
shows the fraction of metals in the cluster at present day as 
a function of redshift they were first incorporated into the 
ICM. The other two panels show the fraction of metals as 
a function of the baryonic mass (middle) and the circular 
velocity (bottom) of the ejecting galaxies. 

The expected delay in the chemical enrichment of the 
ICM is evident if one compares the ejection and wind mod- 
els with the retention model in the top panel of the figure. 
Note that a similar delay is also seen in the metallicity of the 
hot component although it is very small because the reduc- 
tion in total accreted ICM almost balances that in accreted 
metals. Over the redshift range 0-2 the difference between 
the metallicity of the hot component in the retention and 
ejection schemes is only 1.5 per cent solar. As expected, the 
behaviour of the wind model is intermediate between the 
ejection and retention schemes. We find that 60-80 per cent 
of the metals are incorporated into the ICM at redshifts 
larger than 1, 35-60 per cent at redshifts larger than 2 and 
20-45 per cent at redshifts larger than 3 (we obtain lower 
values for the ejection scheme, higher values for the retention 
scheme and intermediate values for the wind model). 

The middle panel shows that 43-52 per cent of the met- 
als are ejected by galaxies with total baryonic mass less than 
1 x 10 10 /i _1 M Q . The distribution of the masses of the eject- 
ing galaxies is approximatively independent of the feedback 
scheme. Although low mass galaxies dominate the luminos- 
ity function in terms of number, they do not dominate the 
contribution in mass. Approximately half of the contribu- 
tion to the chemical pollution of the ICM is from galaxies 
with baryonic masses larger than 1 x 10 10 h~ 1 M@. 

In our model, the star formation and feedback processes 
depend on the virial velocity of the parent substructure. In 
the bottom panel, we show the metal fraction as a function 
of the virial velocity of the ejecting galaxies. We find that 
80-88 per cent of the metals were ejected by galaxies with 
virial velocities less than 250kms _1 . The 'dip' visible for 
the wind model in this panel corresponds to the sharp value 
of 14 r it we adopt for this feedback scheme. Note that the 
results we find from this analysis are very similar to the 
results found bv IKauffmann fc Charlod Jl998l) . even though 
the modelling details are different. 



7 OBSERVATIONAL TESTS FOR DIFFERENT 
FEEDBACK SCHEMES 

Our analysis has shown that properties of galaxies at low 
redshifts are rather insensitive to the adopted feedback 
scheme after our free parameters have been adjusted to ob- 
tain a suitable overall normalisation. 

The observational tests that would clearly distinguish 
between the different models are those that are sensitive to 
the amount of gas or metals that have been ejected outside 
dark matter haloes. 

X-ray observations directly constrain the amount of hot 
gas in massive virialized haloes. The gas fraction tends to 
decrease as the X-ray temperat ure of the system goes down. 
iDavid. Jones fc Formanl lll995l) estimated that the gas-to- 
total mass fraction decreases by a factor ~ 2-3 from rich 
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Figure 9. Fraction of metals, today present in the ICM, as a 
function of the redshift they were first incorporated into the intra- 
cluster gas (upper panel), as a function of the total mass of the 
ejecting galaxy (middle panel) and as a function of the virial 
velocity of the parent substructure (lower panel). Filled circles 
are used for retention model, empty circles are used for ejection 
model and filled squares are used for wind model. 



clusters to groups. In elliptical galaxies, the hot gas frac- 
tion is ten times l ower than in rich cluste rs. These results 
were confirmed by ISanderson et alJ (I2003T) in recent study 
of 66 clusters and groups with X-ray data. One caveat is 
that the hot gas in groups is detected to a much smaller 
fraction of the virial radius than in rich clusters, so it 
is not clear whethe r current estimate s accurately reflect 
global gas fractions llLoewensteinl feoOO ) . Another argument 
for why the gas fraction must decrease in galaxy groups 
comes from constraints from the soft X-ray background. If 
galaxy groups had the same gas fractions as clusters, the 
observed X-ray background at 0.25 keV would exceed the 
observed upper limit by an order of magnitude (|PerJll99sl : 
IWu. Fabian fc NulsenlEoOlh . 

In Fig. I1UI we plot the gas fraction, i.e. Mhot/(Afhot + 
Mvir), as a function of virial mass for all the haloes in 
our simulation. The mass of the hot gas is computed us- 
ing Eq. |5] The results are shown for the case of a short 
re-incorporation timescale (7 = 0.3) in the left column and 
for a long re-incorporation timescale (7 = 0.1) in the right 
column. In this figure we plot results for the M2 simula- 
tion (see Table 0. Note that when we change 7, we also 
have to adjust the feedback efficiency e so that the model 
has the same overall normalisation in terms of the total 
mass in stars formed (for 7 = 0.3, we assume e = 0.2). 
We find that the gas fraction remains approximately con- 
stant for haloes with masses comparable to those of clusters, 
but drops sharply below masses of 10 13 Mq. As expected, 
the drop is more pronounced for the model with a long re- 
incorporation timescale. 

The trends for the wind scheme are similar to those for 
the ejection scheme, but as can be seen, the 'break' towards 
lower gas fractions occurs at lower halo masses, because no 
material is ejected for haloes with circular velocities larger 
than Krit- In the ejection scheme, material is ejected for 
haloes up to a circular velocity of 350 km s . Above this 
value, cooling flows shut down and no stars form in the cen- 
tral galaxies. Galaxies continue, however, to fall into the 
cluster. When a satellite galaxy is accreted, we have assumed 
that its ejected component is re-incorporated into the hot 
ICM. It is this infall of satellites that causes the gas frac- 
tions to saturate in all schemes at halo masses larger than a 
few times 10 13 Af Q . 

The ejection model with a long re-incorporation 
timescale is perhaps closest to satisfying current observa- 
tional constraints. It reproduces the factor 2-3 drop in gas 
fraction between rich clusters and groups. By the time one 
reaches haloes of 1O 12 M0, the gas fractions have decreased 
to a few percent. This may help explain why there has so 
far been a failure to detect any diffuse X-ray emission from 
haloes around late -type spiral galaxies jBenson et aljl2000l : 
iKuntz et al1l2003l) . Note that in the retention model, we 
have assumed that the material reheated by supernovae ex- 
plosions never leaves the halo. This translates into a hot gas 
fraction that is almost constant as a function of M2oo- This 
model is therefore not consistent with the observed decrease 
in baryon fraction from rich clusters to galaxy groups. 

Another possible way of constraining our different feed- 
back schemes is to study what fraction of the metals reside 
in the diffuse intergalactic medium, well away from galax- 
ies and their associated haloes. We now study how metals 
are partitioned among stars, cold gas, hot halo gas and the 
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Figure 10. Fraction of hot gas as a function of M200- Solid circles represent the median and the error bars mark the 20th and 80th 
percentiles of the distributions. 



'ejected component' and we show how this evolves as a func- 
tion of redshift. 

We plot the evolution of the metal content in the dif- 
ferent components for our three different feedback schemes. 
For the ejection scheme we show results for two different 
re-incorporation timescales. The metal mass in each phase 
is normalised to the total mass of metals in the simulation 
at redshift zero. This is simply the yield Y multiplied by 
the total mass in stars formed by z = 0. Note that in the 
retention scheme, all the reheated gas is put into the halo 
and there is no ejected component. 

Fig. llll shows the evolution of the metallicity for an av- 
erage 'field' region of the Universe (our simulation M2). In 
contrast, Fig. 1121 shows results for our cluster. In this plot 
we only consider the metals ejected from galaxies that reside 
within the virial radius of the cluster at the present day. We 
also normalise to the total mass of metals inside this radius, 



rather than the total mass of metals in the whole box. Be- 
cause very little material is ejected from massive clusters in 
any of our feedback schemes, the ejected component always 
falls to zero by z = in Fig. 1121 

Let us first consider the field simulation (Fig. llH . In the 
retention scheme, more than 70 per cent of the metals are 
contained in the hot gas. About 20-30 per cent of the metals 
are locked into stars and around 10 per cent of the metals 
are in cold gas in galaxies. In the ejection scheme, a large 
fraction of the metals reside outside dark matter haloes. This 
is particularly true for the 'slow' re-incorporation scheme, 
where the amounts of metals in stars, in hot gas and in 
the ejected component are almost equal at z = 0. For the 
'fast' re-incorporation model, the metals in the hot gas still 
dominate the budget at low redshifts. 

In the wind scheme, the amount of metals outside viri- 
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alized haloes is considerably lower. There are a factor 2 more 
metals in the hot gas than there are in stars at z — 0. 

At higher redshifts, the relative fraction of metals con- 
tained in cold galactic gas and in the ejected component 
increases. This is because galaxies are less massive and re- 
side in dark matter haloes with lower circular velocities. As 
a result, they have higher gas fractions (see Sec. 14.211 and are 
able to eject metals more efficiently. In the ejection scheme, 
the material that is reheated by supernovae is always ejected 
from the halo, irrespective of the circular velocity of the sys- 
tem. In this scheme, metals in the ejected component (i.e. in 
diffuse intergalactic medium) dominate at redshifts greater 
than 2 in the case of fast re- incorporation and at all redshifts 
in the case of slow re-incorporation. In the wind scheme, ma- 
terial is only ejected if the circular velocity of the halo lies 
below some critical threshold. In this scheme, there are al- 
ways more metals associated with dark matter haloes than 
there are in the ejected component. 

Turning now to the cluster (Fig. I12t . we see that the 
ejected component is negligible in all three schemes. The 
amount of metals outside dark matter haloes never exceeds 
the amount of metals locked up in stars up to redshift ~ 
3. This is because dark matter haloes collapse earlier and 
merge together more rapidly in the overdense regions of the 
Universe that are destined to form a rich cluster. Although 
metals may be ejected, they are quickly re-incorporated as 
the next level of the hierarchy collapses. Note that the metals 
in the cold gas also fall sharply to zero at low redshifts. 
This is because galaxies that are accreted onto the cluster 
lose their supply of new gas. Stars continue to form and the 
cluster galaxies simply run out of cold gas. 

Comparison of Fig. llll with Fig. ll2l suggests that cosmic 
variance effects will turn out to be important when trying to 
constrain different feedback schemes using estimates of the 
metallicity of the intergalactic medium deduced from, for 
example, CIV absorptio n systems in the spectra of quasars 
Je.g. Schave et al.l2003h . Nevertheless, the strong differences 
between our different feedback schemes suggest that these 
kinds of measurements will eventually tell us a great deal 
about how galaxies ejected their metals over the history of 
the Universe. 

Finally, in Fig. ll3l we show the evolution of the different 
phases for the slow ejection scheme, that can be considered 
as our 'favourite' model. In the left panel we plot the evo- 
lution in mass of the different phases for all the galaxies 
within the virial radius in our cluster simulation. In agree- 
ment with what is shown in Fig. |S] the stellar component 
grows very slowly after redshift ~ 3, when most of the stars 
in the cluster have already been formed. The cold gas and 
the ejected mass decrease as galaxies are accreted onto the 
cluster and stripped of their supply of new gas. The hot gas 
mass increases because of the accretion of 'diffuse material'. 
In the right panel we plot the evolution of the mass density 
in each reservoir for our field simulation. Note that here the 
evolution in the stellar mass density is more rapid than in 
the field. 



8 DISCUSSION AND CONCLUSIONS 

We have presented a semi-analytic model that follows the 
formation, evolution and chemical enrichment of galaxies in 
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Figure 11. Evolution of the metal content of different phases for 
the three different models used in this paper and for a typical 
'field' region. In each panel the solid line represents the evolu- 
tion of the metal content in the stars, the dashed line the cold 
gas, the dashed-dotted line the hot gas and the long-dashed line 
the ejected component (not present in the retention model). The 
metal content in each phase is normalised to the total mass of 
metals produced from all the galaxies considered. 
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Figure 13. Evolution of the different phases for the slow ejection scheme. The left panel shows the results for the galaxies within the 
virial radius of our cluster simulation, while the right panel shows the results for a typical 'field' region. 



a hierarchical merger model. Galaxies in our model are not 
closed boxes. They eject metals and we track the exchange 
of metals between the stars, the cold galactic gas, the hot 
halo gas and an ejected component, which we identify as the 
diffuse intergalactic medium (IGM). 

We have explored three different schemes for imple- 
menting feedback processes in our models: 

• in the retention model, we assume that material re- 
heated by supernovae explosions is ejected into the hot halo 
gas. 

• In the ejection model, we assume that this material is 
ejected outside the halo. It is later re-incorporated after a 
time that is of order the dynamical time of the halo. 

• In the wind model, we assume that galaxies eject ma- 
terial until they reside in haloes with V v i r > Krit- Ejected 
material is also re-incorporated on the dynamical time-scale 
of the halo. The amount of gas that is ejected is proportional 
to the mass of stars formed. 

In all cases we can adjust the parameters to obtain rea- 
sonably good agreement between our model predictions and 
observational results at low redshift. The wind scheme is per- 
haps the most successful, allowing us to obtain remarkably 
good agreement with both the cluster luminosity function 
and the slope and zero-point of the Tully-Fisher relation. 
All our models reproduce a metallicity mass relation that is 
in striking agreement with the latest observational results 
from the SDSS. By construction, we also reproduce the ob- 
served trend of increasing gas fraction for smaller galaxies. 
The good agreement between models and the observational 
results suggests that we are doing a reasonable job of track- 
ing the circulation of metals between the different baryonic 
components of the cluster. 

We find that the chemical enrichment of the ICM oc- 
curs at high redshift: 60-80 per cent of the metals are ejected 
into the ICM at redshifts larger than 1, 35-60 per cent at 
redshifts larger than 2 and 20-45 per cent at redshifts larger 



than 3. About half of the metals are ejected by galaxies with 
baryonic masses less than 1 x 10 10 h~ 1 MQ. The predicted 
distribution of the masses of the ejecting galaxies is very 
similar for all 3 feedback schemes. Although small galaxies 
dominate the luminosity function in terms of number, they 
do not represent the dominant contribution to the total stel- 
lar mass in the cluster. Approximately the same contribution 
to the chemical pollution of the ICM is from galaxies with 
baryonic masses larger than 1 x 10 10 h~ 1 M.Q. 

Finally, we show that although most observations at 
redshift zero do not strongly distinguish between our dif- 
ferent feedback schemes, the observed dependence of the 
baryon fraction on halo virial mass does place strong con- 
straints on exactly how galaxies ejected their metals. Our 
results suggest that gas and its associated metals must be 
ejected very efficiently from galaxies and their associated 
dark matter haloes. Once the material leaves the halo, it 
must remain in the diffuse intergalactic medium for a time 
that is comparable to the age of the Universe. 

Future studies of the evolution of the metallicity of the 
intergalactic medium should also be able to clarify the mech- 
anisms by which such wind material is mixed into the envi- 
ronment of galaxies. 
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